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ABSTRACT 


The use of an analog computer in the simulation of a dc arc and its control system 
is described. The detailed model design is presented. Control parameters determined 
with the model are shown to be optimum with the real system. Repeatability of system 
control to a programmed input to ±1 percent is achieved for both the model and the real 
system. Data from the model and the arc system show close agreement except in some 
areas of arc instability. 
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SUMMARY 

The use of an analog computer in the simulation of a dc arc and its control system 
is described. The detailed model design is presented. Control parameters determined 
with the model are shown to be optimum with the real system. Repeatability of system 
control to a programmed input to ±1 percent is achieved for both the model and the real 
system. Data from the model and the arc system show close agreement, except in some 
areas of arc instability. 


INTRODUCTION 

A new technique (ref. 1) in emission spectrochemical analysis being developed at 
the NASA Lewis Research Center requires control of the emission output of the static 
argon stabilized dc arc. Previously, constant- current control of the arc was used in 
spectrochemical analysis to improve the reproducibility of the sample vaporization, and, 
hence, the analytical precision. For further improvement in the precision of quantita- 
tive analysis, the use of a current control system based on the arc light intensity was 
deemed necessary. A control system (ref. 1) was incorporated to control the dc arc 
current such that the dispersed emission intensity received at a chosen photomultiplier 
(i. e. , a selected element emission line) would follow a given programmed function of 
time repeatably to within 1 percent. This programming of the arc emission results in 
an indirect control of the arc temperature and plasma characteristics. Since the param- 
eter of primary importance was the integral of the light intensity over the program time, 
transient errors of long duration could not be tolerated. However, difficulties were in- 
curred in attempting to achieve stable operation of the control loop at the required gains. 
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The nonlinear and time-varying characteristics of the arc and its control system made 
theoretical analysis of the system very complex. 

The use of analog computer simulation techniques appeared suitable and advanta- 
geous to the analysis of the control system. Nonlinear and time-varying parameters 
could be incorporated readily in the computer model. If successful simulation could be 
attained, the computer model might be useful in determining the control parameters for 
the system as a function of the electrodes to be used, the spectral line to be controlled, 
and the programmed input to be followed. Both time and cost might be saved over trial- 
and-error stabilization methods using numerous samples in the dc arc. Increased un- 
derstanding of the arc control system and of the arc processes themselves might pos- 
sibly be obtained by study of the simulator. Even without perfect simulation, areas 
could be singled out for further study by the elimination from consideration of the well- 
described characteristics of the arc system. 

A study of the literature to determine if any analog simulation of dc arc emission 
had been reported was unfruitful. A literature search for mathematical models of the 
arc emission process was also unsuccessful. Studies have been made of the thermal and 
radiation characteristics of a static argon stabilized dc arc (ref. 2); but these studies 
are not directly applicable to spectral line emission. The emission intensity of a given 
element line is a complex process sensitive to temperature changes in the arc and, 
therefore, difficult to predict. The arc model used in this report was, therefore, based 
solely on experimental data obtained with the specific static argon dc arc used in the 
spectrochemical analysis. 

This report describes the computer model used to simulate the arc and its control 
system and the information derived from the simulation. The use of simulation tech- 
niques required the experimental determination of the system functions. Approximations 
using experimental data were often required to permit a reasonably concise, analytical 
description of the system characteristics. With this knowledge of the arc and its control 
system mathematically stated, an analog computer model of the system was developed. 
The computer methods used and the results obtained are discussed in this report. Com- 
parisons between the model and the real system were used to evaluate the simulation. 
Repeatability of control to a programmed input to ±1 percent was achieved. Parameters 
determining stability of the model were employed with the real system, noise effects 
were studied, and real system changes were incorporated based on the computer model 
results. 


ARC CONTROL SYSTEM 

The arc used in emission spectrochemical analysis at the Lewis Research Center is 
a static argon stabilized dc arc (ref. 3). The major components of the arc are shown in 
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Figure 1. - Static argon stabilized dc arc. 


Figure 2. - Block diagram of dc arc and control system. 


figure 1. The arc anode is fabricated with a cylindrical cup for holding the sample to be 
vaporized. The cathode is tipped with a tantalum bead. The emission from the arc gap 
is optically coupled to the spectrometer for analysis. The arc is contained in a chamber 
for control of the pressure and constituents of the arc atmosphere. 

The control system (ref. 1) used with the dc arc achieves control of the arc emis- 
sion intensity output by regulating the current to the arc. The controlled parameter is a 
selected spectral atom line of silver, Ag (350.2 nm). The block diagram of the arc con- 
trol system is shown in figure 2. The intensity I, as detected at the Ag sensitive photo- 
multiplier in the emission spectrometer, is proportional to the intensity of emission of 
the Ag in the arc gap. The electrometer amplifier converts the current output of the 
photomultiplier to a voltage output proportional to the Ag emission intensity. The 
curve follower provides a time-varying voltage V 2 as a program for emission intensity 
output to which Vj is to be controlled. (The method for obtaining Vg will be described 
in a later section. ) Comparison of with the curve-follower voltage Vg results in 
an error voltage e = Vg - Vj. The control unit operates on e, forming K-^e, K 2 (de/dt), 

and Kg J e dt, whose summation is the control voltage Vg. The control voltage Vg 


controls the thyratron current supply output i by phase shifting the thyratron grid sig- 
nal. A positive change in Vg increases the arc current, while a negative change de- 
creases the current. Current changes in the arc are such as to minimize e and thus 
cause V x to follow the desired program V 2 . A constant voltage V 4 is used to provide 
a minimum arc current independent of e . 

The optimum values for K^, Kg, and Kg to achieve high gain control and stability 
are dependent on the arc function, the light measuring circuitry, and the current supply 

circuitry. The integral control voltage Kg J e dt is used to maintain a steady- state 

error of zero. A control voltage K(de/dt) is used to provide system damping. To 
lessen the effects of high-frequency noise on the system response, the derivative signal 
is low-pass filtered. 
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DETERMINATION OF SYSTEM EQUATIONS 


Analog computer simulation of the dc arc and the control system requires formula- 
tion of a mathematical model describing the system in terms of the function of its com- 
ponents. 

The block diagram of the system (fig. 2) will be described by three functions exper- 
imentally determined from measurements on the arc and the associated electronics. A 
function Vj(i, t) will be used to describe the voltage output corresponding to the inten- 
sity of the silver emission, with the arc current and time as the independent variables. 
The function Vg(e, t) is the functional operation of the control unit acting on the error 
signal. Finally, i(Vg, t) corresponds to the arc current as a function of the control unit 
output voltage and time. 


Function V]0,t) 

Description of Vj(i, t) required experimental data of the steady-state and dynamic 
response of the Ag spectral atom line intensity to the arc current. 

These experimental data showed that could be described by the equation 

V l (i,t)=f 13 {f u [f 11 (i,t)]} (1) 

where this expression uses functional notation (i. e. , f ^ is a function of f jj, and f^g 
is a function of f ^). The function relating the arc emission intensity to the arc current 
is f jj. The function f ^ describes the transient response of the emission intensity to 
changes in arc current, and f^g is the transient response of the electrometer amplifier. 

The function f^j was determined by measuring the Ag emission intensity in terms 
of Vj for several arc current values. The time dependence of f ^ is a result of the 
arc temperature initially increasing at ignition and also of the decreasing amount of 
sample. For a constant initial sample size, therefore, f ^ is a current- dependent, 
time -varying function. The variation in silver emission intensity from sample to sample 
(as large as 20 percent) made it necessary to average the results of several samples at 
each current level. 

The peak light output during the vaporization of the sample was plotted as a function 
of arc current squared. Figure 3 shows this plot for two types of arc anode. In the re- 
gion of interest (20 to 38 A of arc current) the data can be approximated by a linear 
function of the variable current squared. The peak Ag emission intensity for the anode 
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Figure 3. - Peak silver emission intensity as function of arc current. 


TABLE I. - VALUES OF CONTROL 


SYSTEM CONSTANTS 


Constant 

Value 

A 

7xl0 -6 V/A 2 

B 

3. 6X10 2 A 2 

C 

900 A 2 

D 

17 A 

E 

38 A 

F 

30 A/V 

T 1 

1. 4 sec 

r 2 

0. 2 sec 

t 3 

0. 1 sec 


of primary interest in spectrochemical analysis (type M) is described by the function 

V lp = A(i 2 -B) (2) 

(Table I lists the values of the constants used in this study. ) 

The Ag emission intensity output for 30 amperes constant current, with a 4- 
milligram sample of silver chloride (AgCl), is shown as a function of time in figure 4. If 
the curves obtained for other current levels are normalized along their time axes to the 
30-ampere curve on the basis of the ratio of the current squared to 900 amperes squared, 
all of the normalized curves approximate in shape the 30-ampere Ag curve. (This is a 
reasonable result, since the integral of the Ag line intensity over the entire sample va- 
porization should be approximately constant for a constant sample size. ) Therefore, the 
time function of the Ag light output can be described by 

f U= V l,p f < t '> (3) 

where 
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Figure 4. - Approximate silver emission intensity for 4-milligram sample of 
silver chloride at arc current of 30 amperes. (Test conditions from ref. 3. ) 


C = 900 amperes squared 


and V, f (t T ) is as shown in figure 4 for 


i 2 


900. Therefore, 


2 /i 2 \ 

fii(i,t) = A(i 2 - B)f ( — 1 1 

\ c / 


(4) 


To obtain the transient response of the Ag light to current, step changes in current 
were applied to the arc during its operation and the Ag light response was measured. 

The Ag line was found to approximate a first-order (single time constant) response with 
a time constant Tj (where was determined to be 1. 4 sec for electrode type M). To 
within the ±20 percent error in measurements, was found to be independent of i. 
Therefore, f^ can be described by the differential equation 


af 


12 


at 


+ f 12 " f ll 


(5) 


The electrometer amplifier was known to contribute another first-order term with a 
time constant T 2 (0.2 sec). Since the amplifier function is linear, its time response 
f can also be described by an independent first-order differential equation 


df ig 

r 2 — + f 13 = f 12 


( 6 ) 
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The overall function Vj(i, t) is obtained from equations (1), (4), (5), and (6) and can 
be described by the second-order differential equation 


t 1 t 2 


3 2 V< dV, 

-i+( T l + T 2 ) + Vj 

at 2 1 2 3t 


A(i 2 - B) f 



( 7 ) 


Function V 3 

The function Vg(e, t) for the control unit is related to the sum of the amplifier out- 
puts generating the proportional, integral, and derivative signals by the equation 


V 3 (t) = Ki€ + |£ + Kj £ € 


dt 


( 8 ) 


The values for the gains (K^, Kg, and Kg) are separately adjustable. 


Function i 

The phase shifting and thyratron circuitry of the arc current supply are described 
by the function i(Vg,t). A bias voltage V4 (see fig. 2) is provided to maintain a mini- 
mum current level D of about 17 amperes. The maximum current E that can be sup- 
plied through the thyratron circuits is 38 amperes. The arc current supplied as a func- 
tion of voltage Vg is a nonlinear function for currents near the 38-ampere upper limit. 
However, the ampere-per-volt gain can be closely approximated by a linear relation and 



Figure 5. - Linear approximation of arc current as 
function of control voltage. 
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a limit at 38 amperes, as shown in figure 5. The linear slope F is 30 amperes per 
volt. With consideration of the bias current, the function is described by 


i = D + FVr 


[for V, 


E - 


5 ) 


(9a) 


i = E 


'for Vo > 
t o 


E - D 


(9b) 


The load inductance filtering the thyratron current waveform results in a first-order re- 
sponse of the current to the control voltage. The time constant Tg = 0. 1 second. 

The function i(Vg, t) may now be represented by the equations 


Tg — + i = D + FVg 
3t 


(■ 


for V 3 < 


E - D' 


(10a) 


i = E 


(*“ V 3>^) 


(10b) 


This completes the mathematical model of the control system necessary for simulation 
techniques. 


Effect of System Changes on Model 

Although this simulation of the arc control system concentrates on the use of an Ag 
spectral line as its control parameter, it is entirely possible that other lines may be 
used in other applications. It was experimentally determined that the equation describ- 
ing the arc emission does not have the same form for different element emission lines. 

Preliminary data on the emission of a potassium K (404. 7 nm) line showed that the 
arc function was of the form 


V x (t) = gl (i, t) = g 14 ^g 13 | g 12 [g n (i, t)]}) 

where g 14 (i, t) describes the natural profile of the K emission, g 14 involves the am- 
plifier time constant, and and gjg are two functions involving the transient re- 
sponse of the light emission process and are described by 


8 



3 ®12 

T l~^- + ^12 ~ g 13 

The response to a step change in arc current is a rapid increase to an output greater 
than the steady-state output followed by a first-order decay to the steady- state output. 

Therefore, simulation of the system performance based on the control of other 
spectral lines would require significant changes in the system description. Experi- 
mental data on the response of a given spectral line to arc current at constant sample 
size must be obtained. 

Similarly, as was illustrated by figure 3, the anode used in the dc arc affects the 
spectral emission of silver. Anode N achieved higher Ag emission intensity for a given 
current than did anode M. Though the transient response of the arc to current changes 
with anode N also exhibited a first-order response, the time constant was measured as 
1. 8 seconds, slightly greater than that for anode M. 

These examples of changes in the system illustrate the need for experimental data 
to determine a model for the complex emission processes of the dc arc. 


ANALOG SIMULATION OF DC ARC AND CONTROL SYSTEM 

With a sufficient mathematical description of the arc and the control system com- 
pleted, the simulation of the system by analog computer follows directly from the sys- 
tem equations (7), (8), (10), and 


e = V 2 - V t (11) 

Equation (7) may be realized on an analog computer by the familiar method of solving for 

2 9 

the highest derivative 3 V-j/St and using computer integrators to form the lower order 
terms. 



(t i + To) 3V i -i i o / «2 \ 

_1 JL. i L_ Vj + l - A(i 2 - B)f(i-t) 

t 1 t 2 9t t 1 t 2 t 1 t 2 L ' C ' 


( 12 ) 


Equation (8), involving the variable i, is similarly solved for the highest order deriva- 
tive. 
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(13) 


^i=-J- + -L(D + FVo) 
at r 3 t 3 

o 

The functions V 2 (t) and f(i t/C) are readily formed on an analog computer with 
diode function generators. The function generators present an output voltage as a func- 
tion of an input voltage. If this input voltage is linearly increasing with time, the output 
voltage will be a function of time. To generate f(i^t/C), the input voltage to the func- 
tion generator is the time integral of the term bi , where b is a constant adjusted such 
that a current of 30 amperes results in the function f(t) in real time. An increase in i^ 
will proportionally increase the voltage ramp to the function generator, and, hence, de- 
crease the time required to reach a given point in the output function. 

The Ag emission intensity function which the output Vj(t) is to follow is V 2 (t); this 
function is shown in figure 6. The steady- state level is the Ag emission intensity outfxit 
generated by an arc current of 30 amperes. The function V 2 was chosen to control the 
Ag emission to values close to the natural constant current emission at 30 amperes. 

This restraint on V 2 was based on the advisability of reproducing arc emission charac- 
teristics for which previous experience in spectrochemical analysis had been obtained. 
The choice of V 2 also guaranteed that the analytical model of the arc would best fit the 
operating conditions with the control system, since the model was derived from experi- 
mental data at constant arc current. A deviation from the natural emission profile was 
programmed for the start of Ag emission. 

The ramp portion of the curve precedes in time the natural profile of the Ag output 
at ignition of the arc at a current level of 30 amperes. This will result in a large cur- 



0 4 8 12 16 20 24 28 

Time, sec 

Figure 6. - Program for silver emission intensity output. 
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rent signal at the initial portion of the ramp to cause the earliest possible emission of 
the Ag at ignition. This criterion is conducive to the requirements for spectrochemical 
measurements of other volatile elements in the sample. It is also desirable, however, 
to limit the initial current surge below the maximum thyratron current, so as to main- 
tain the current in control throughout the vaporization of the entire sample. 

The slight nonlinearities in the shape of VgOO at maximum voltage (as shown in 
fig. 6) are a result of limitations in the use of diode function generators for connecting 
linear segments of a function. The effects on the model response are negligible. 

The analog computer simulation of the control system was achieved on a 10-volt 
full-scale analog computer. Since the time parameters involved were suitable for real 
time programming, no time scaling was necessary. However, voltage amplitude scaling 
was necessary in order to observe the variables of interest: Vj(t), V^t), e(t), and i(t). 
Table H lists the computer variables and magnitudes, and the computer scale factors. 


TABLE n. - COMPUTER VARIABLES 
AND UNITS 


Variable 

Maximum value 

Scale factor 

V^t) 

10 mV 

1 V/mV 

V 2 (t) 

10 mV 

1 V/mV 

iW 

38 A 

0. 1 V/A 

i 2 (t) 

1444 A 2 

0. 001 V/A 2 


Rewriting equation (13), and expressions for i(t) and Vg(t), for the computer vari- 
ables results in the computer equations 



V (t) = 5k (10 3 e) + (lO 3 h\ + J* 
10 3 10 3 ' 9t ' lo 3 0 


10 3 e dt 


(16) 
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Figure 7. - Analog computer model for dc arc and control system. 
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Figure 7 shows the complete computer simulation of equations (14), (15), and (16). The 
computer variables are shown in the diagram where they occur, although the amplitude 
scaling is not explicitly labeled. The function generators are labeled as to the function 
programmed. The symbols used for integrator, summer, operational amplifier (no in- 
ternal feedback), and function generators are listed in the key. The numbers on the input 
lines to the integrators and summers represent the gain for the respective inputs (i. e. , 
gain of 10 or 1). 

The differentiating circuit provides an approximate first derivative of the error sig- 
nal by 


0V 

v «rt< t > = K af- - a -^r < 17 > 

where a is a variable set at 0. 02 to limit the derivative signal at frequencies above 
8 hertz, thereby lessening the effect of high-frequency noise on the circuit response. 

The potentiometer settings are shown in the diagram (fig. 7) as percent of signal 
transmitted, with the gain potentiometers labeled as a function of proportional, integral, 
or derivative gain. 

Provision for the simulation of arc noise is provided at the inputs labeled ” Pulse.” 
The natural profile of the Ag emission line as a function of time, shown in figure 4, is 
an average curve which does not show the higher frequency noise pulses encountered in 
the arc. This noise is simulated by the manual pulses provided by external circuitry. 
The pulse rise time is in the millisecond range, while the decay time constant is 
0. 1 second. 


RESULTS AND DISCUSSION 

The analog computer simulator was used to determine the optimum proportional, 
integral, and derivative gains for the real system. Study of the simulator operation re- 
sulted in improvement of the real system design. Comparisons of simulator and real 
system responses provided information on the accuracy of the arc model. 

Gain parameters exhibiting stable control on the model were used with the real sys- 
tem and comparisons were made by using the real system and the model responses to an 
emission intensity program. The bases of comparison are the system stability, the 
emission-intensity waveform, and the arc- current waveform. A strict comparison in 
waveforms is limited by the sample-to-sample variations in the Ag emission intensity 
for equal quantities of AgCl and by limitations in accurately simulating the random fluc- 
tuations in emission intensity. It was noted earlier that the magnitude of the emission 
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intensity output for constant arc current was found to vary from sample to sample by as 
much as 20 percent. 


Simulator Response 

The simulator characteristics were studied for various magnitudes and interrela- 
tions of proportional, integral, and derivative gains. The proportional gain was in- 
creased in steps until only minor improvements in control were obtained or until unsta- 
ble control occurred. The integral gain was adjusted to reduce the initial error as rap- 
idly as possible without excessive overshoot of the emission intensity program and with- 
out causing the arc current to reach the current limit. 

The effect of integral gain is illustrated by figure 8. The integral gain Kg for 
curve 2 is twice that for curve 1 (1000:500). The emission intensity output V^, the pro- 
grammed input Vg, the arc current i, the natural profile function f(t'), and the error 
voltage e are displayed as recorded in real time. The differences in i and e for 
curves 1 and 2 are to be noted. The increase in Kg for curve 2 results in a higher ini- 
tial surge in i at t = 5 seconds, resulting in a decrease in e to zero at about 
t = 6 seconds and a subsequent overshoot (negative error signal) until t = 9 seconds. 

For curve 1, the error does not decrease to zero as rapidly, but no significant overshoot 
occurs. (Although the computer runs are continued to about t = 24 sec, sample vapor- 
ization in the real system would be terminated at about 22 sec; this is the time at which 
Vg starts to decrease. ) 

The derivative gain was adjusted to provide system damping for minimization of 
overshoot and oscillation in the system waveforms for ideal conditions (as in fig. 8) and 
also for noise disturbances applied to the system. 

Figure 8 shows a slight anomaly in the waveforms f(t T ) and Vj. During the time 
t < 4 seconds, Vj remains zero while f(t’ ) slowly increases because of the minimum 
bias current of 17 amperes. The voltage is zero because the bias current is insuf- 
ficient to result in a positive value for V-, as determined by the linear approximation 

*9 P 

in figure 3. The real dc arc does exhibit a slight Ag emission for minimum arc current, 
but the integral of this value of emission for a 4-second period is negligible compared to 
the total Ag emission. 

Since arc noise disturbances in the form of spikes in emission intensity were always 
present in the real system operation, the computer studies were concentrated on deter- 
mining the effects of noise as a function of the control gains. Noise pulses were applied 
at f(t') and at (see fig. 7). The pulses at f(t') were to simulate discontinuities in 
the vaporization of the AgCl at the arc anode. Inconsistencies in the form or purity of 
the AgCl could conceivably cause this effect. Pulses were applied at to simulate 


14 












40 

30 

20 

10 

0 





Figure 11. - Model response to noise pulses 
applied at Vp with high derivative gain. 
Proportional gain, Kj, 1000; derivative 
gain, K 2 , 1500; integral gain, K 3 , 500. 




system operation under noise conditions representing changes in the arc chamber or in 
the spectrometer and electronics which would affect the Ag emission intensity output di- 
rectly. These changes might be the result of arc temperature fluctuations, minor plas- 
ma displacements, or electronic noise. 

Figure 9 shows the operation of the simulated arc system with noise disturbances at 
f(t’), with optimum integral and derivative gains for the proportional gain of 1000. After 
the initial transient error, control to within ±1 percent of program is achieved. The 
transient effects of the current changes required to correct for the noise disturbances 
are negligible. 

Figure 10 illustrates the system operation for disturbances at Vj. These disturb- 
ances are more disruptive to system control than were pulses at f(t’). The control 
gains are adjusted for minimal transient effects - the transients reduce to near zero 
within 1 second after the disturbance. The effect of a large change in derivative gain 
(1500) for disturbances at Vj is shown in figure 11. The current oscillations result in 
large (10 percent) errors of several seconds duration. 


Comparisons of Model and Real System 

The use of optimum control gain parameters, derived from the model study, gener- 
ally resulted in stable, optimum control with the real system. For proportional and in- 
tegral gains of 1000 and 500, respectively, stable arc control was achieved with a deriv- 
ative gain of 300. The real system response under these gain conditions is illustrated 
by figure 12. Control to the program to within 1 percent is achieved, with the exception 
of the repeatable initial lag. The current waveform is similar to that achieved on the 
simulator (as in fig. 8). The emission intensity disturbance halfway up the ramp portion 
of the program is corrected by a current change whose transient decays within about 
2 seconds. (In fig. 12, because of the recording pen offsets, the current trace leads the 
emission intensity output trace by 1 sec. For that reason, the output disturbance does 
not exactly coincide in time with the current response as recorded in the figure. ) 

Although simulator response did not vary appreciably (with Kj = 1000 and Kg = 500) 
for a range of derivative gains of ±30 to 40 percent, this was not the case with the real 
system. Variation of the derivative gain of more than 20 percent from optimum resulted 
in unstable operation of the control system. It appears that instabilities in the dc arc 
are triggered by current oscillations which result from noise disturbances. This dis- 
crepancy between the model and real system responses will be discussed later. 

Comparisons of the arc system and its model were also made at higher proportional 
and integral gains. In general, the derivative gains for optimum damping were about 
equal. For proportional gains of 5000 on the model, the band of derivative gains for 
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Figure 12. - Real system response to 
input program. 


stable response to noise disturbances was about 20 percent, and even the best responses 
exhibited transient oscillations, as shown in figure 13. It can be noted in the figure that 
although the initial transient error for t < 4 seconds has been reduced, the high control 
gains have caused the initial current surge to reach the current limit at 38 amperes. 
During this time of about 1 second duration, the current control was saturated and could 
not respond to noise disturbances. Although stable control of the dc arc could also be 
attained at this control gain, no significant improvement in Ag emission control was 
achieved because of more pronounced noise transient effects. Stable operation was not 
achieved at proportional gains of 10 000 or above for either the model or the real system. 

Some improvement in the dc arc control system resulted from studies on the com- 
« puter model. Measurements of the derivative control gain for the system prior to simu- 

lation showed that the derivative gain decreased as a first-order function above a fre- 
quency of 3 hertz. Study of the model responses showed that improvements in system 
stability could be obtained if the derivative gain corner frequency was extended to 
8 hertz. A change was made in the real system to incorporate this characteristic. The 
result was a significant improvement in real system stability. 
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Real System Instabilities 


The analog simulator exhibited limitations in accurately reproducing arc instabili- 
ties which can be disruptive to control of the emission intensity. 

The foremost example of the model limitations involved arc instabilities related to 
the alinement of the arc electrodes. This little-understood phenomenon became evident 
during several of the initial sample vaporizations with the real arc system. With gain 
and arc operating conditions for which stable arc operation had been obtained regularly, 
control of the Ag emission line became unstable. These instabilities at first seemed to 
be related to the cathodes used (ref. 3), since insertion of a new cathode often resulted 
in control oscillations. However, it was then determined that reinsertion of previously 
stable cathodes did not necessarily guarantee stable control. It was determined experi- 
mentally in open-loop operation that there were no apparent changes in the arc function 
f^ for stable or unstable cathodes. 

An explanation of this phenomenon was arrived at through subsequent testing. It 
was determined that a slight misalinement of the axes of the arc cathode and anode re- 
sulted in stable operation. Since the cathode construction techniques (ref. 3) often re- 
sulted in slightly unsymmetrical cathodes, reinsertion of a cathode did not necessarily 
place it in the same position relative to the anode. When care was taken to insert the 
cathode in a predetermined orientation, stable control of the arc was achieved. 

A possible explanation of the dependence of arc stability on electrode positioning 
followed from visual observation of the plasma streaming between the arc electrodes. 
For the slight misalinement condition of the electrodes, the plasma streaming from 
cathode to anode appeared to achieve a stable, circular flow pattern. With the cathode 
tip closely alined with the anode axis, however, the plasma flow pattern impinged on the 
anode directly, the pattern was irregular, and arc instabilities occurred. 

It was noted earlier that the range of derivative gains for stable control of the real 
system was narrower than the range of optimum gains for the model. For proportional 
gains of about 5000, large changes in derivative gain from the optimum value were nec- 
essary to cause instability of the model. The real arc appeared to be more susceptible 
to current oscillations caused by noise disturbances. Once unstable, the emission in- 
tensity oscillations were of larger magnitude than would be expected on the basis of the 
computer model studies. 

Further studies of the arc function under various conditions of plasma flow and for 
external and internal disturbances are necessary to provide sufficient information to 
better characterize the arc emission processes for computer simulation. 
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General Observations 


The computer simulation of the dc arc and the control system were used in this 
study as a tool to determine optimum gain parameters for control of the Ag emission in- 
tensity. It has been noted earlier that changes in the element emission line being con- 
trolled or in the arc electrodes necessarily change the analytical description of the sys- 
tem, Further studies with simulation techniques might be helpful in determining what 
characteristics of the parts of the system result in optimum control. Also, the use of 
various input programs could readily be investigated with the model for optimization of 
control with consideration for the requirements of spectrochemical analysis. These 
uses for the analog simulator would all be enhanced by improvements in the analytical 
description of the arc emission processes. 


SUMMARY OF RESULTS 

An analog simulation of a dc arc and its control system was developed and used to 
determine optimum control gain parameters for the real system. The analytical model 
for the arc emission process was derived from experimental data. The following are 
the results from the study of the analog model and from comparisons between the model 
and real system responses: 

1. The model was used to obtain control gain parameters for stable operation of the 
real dc arc system. 

2. The effects of noise disturbances on the control of the model emission intensity 
were studied. 

3. The analytical model of the dc arc was shown to be a function of the dc arc pa- 
rameters, including the controlled element emission line and the arc anode. 

4. The model of the dc arc emission process was found to have some limitations in 
simulating instabilities in the arc. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, April 23, 1968, 

129-03-14-04-22. 
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